Electron spin evolution induced by interaction with nuclei in a quantum dot 
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We study the decoherence of a single electron spin in an isolated quantum dot induced by hyperfine 
interaction with nuclei for times smaller than the nuclear spin relaxation time. The decay is caused 
by the spatial variation of the electron envelope wave function within the dot, leading to a non- 
uniform hyperfine coupling A. We show that the usual treatment of the problem based on the 
Markovian approximation is impossible because the correlation time for the nuclear magnetic field 
seen by the electron spin is itself determined by the flip-flop processes. The decay of the electron spin 
correlation function is not exponential but rather power (inverse logarithm) law-like. For polarized 
nuclei we find an exact solution and show that the precession amplitude and the decay behavior can 
be tuned by the magnetic field. The decay time is given by hN/A, where iV is the number of nuclei 
inside the dot. The amplitude of precession, reached as a result of the decay, is finite. We show that 
there is a striking difference between the decoherence time for a single dot and the dephasing time 
for an ensemble of dots. 

PACS numbers: PACS numbers: 85.35.Be; 73.21. La; 76.20.+q; 76.60.Es 



I. INTRODUCTION 

The spin dynamics of electrons in semiconducting 
nanostructures has become of central interest in recent 
years |l], p|. The controlled manipulation of spin, and in 
particular of its phase, is the primary prerequisite needed 
for novel applications in conventional computer hardware 
as well as in quantum information processing. It is thus 
desirable to understand the mechanisms which limit the 
spin phase coherence of electrons, in particular in GaAs 
semiconductors, which have been shown || to exhibit 
unusually long spin decoherence times Ti exceeding 100 
ns. Since in GaAs each nucleus carries spin, the hyper- 
fine interaction between electron and nuclear spins is un- 
avoidable, and it is therefore important to underdffect on 
the electron spin dynamics Q . This is particularly so for 
electrons which are confined to a closed system such as 
a quantum dot with a spin 1/2 ground state, since, be- 
sides fundamental interest, these systems are promising 
candidates for scalable spin qubits ||. For recent work 
on spin relaxation (characterized by Ti times) in GaAs 
nanostructures we refer to Refs.|| 0, |8|. 

Motivated by this we investigate in the following the 
spin dynamics of a single electron confined to a quan- 
tum dot in the presence of nuclear spins (see Fig. 1)0. 
We treat the case of unpolarized nuclei perturbatively, 
while for the fully polarized case we present an exact 
solution for the spin dynamics and show that the de- 
cay is non-exponential and can be strongly influenced 
by external magnetic fields. We use the term "decoher- 
ence" to describe the case with a single dot, and the term 
"dephasing" for an ensemble of dots [Eof. The typical 



fluctuating nuclear magnetic field seen by the electron 
spin via the hyperfine interaction is of the order of |ll|| 
~ Af^y/Ngns), with an associated electron precession 
frequency ujn ~ A/y/W, where A is a hyperfine constant, 
g the electron g-factor, the Bohr magneton, and N 
the number of nuclei inside the dot. For a typical dot size 
the electron wave function covers approximately N = 10 5 
nuclei, then this field is of the order of 100 Gauss in a 
GaAs quantum dot. The nuclei in turn interact with each 
other via dipolar interaction, which does not conserve the 
total nuclear spin and thus leads to a change of a given 
nuclear spin configuration within the time T„2 ~ 10 _4 s, 
which is just the period of precession of a nuclear spin in 
the local magnetic field generated by its neighbours. 

We note that there are two different regimes of in- 
terest, depending on the parameter u>nt c , where r c is 
the correlation time of the nuclear magnetic field H^v (t) 
seen by the electron spin via the hyperfine interaction. 
The simplest case is the perturbative regime u>nt c <C 1, 
characterized by dynamical narrowing: different random 
nuclear configurations change quickly in time, and, as 
a result, the spin dynamics is diffusive with a dephas- 
ing time ~ l/(u;jyT c ). This case is realized, for example, 
for a system of quantum dots (or shallow donors) when 
the hopping rate 1 /r c of the electron between neighbour- 
ing dots is high. The problem of electron spin relax- 
ation for the case of electron hopping between shallow 
donors in GaAs was studied in Ref. [l^. Using a pertur- 
bative approach, we easily obtain the following formulas 
for the longitudinal (Ti) and transverse (T2) spin relax- 
ation times jd| 
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where lo z — gfiBB/h is the Larmor spin precession 
frequency in the external magnetic field B directed along 
the z-axis, and < ... > means the ensemble average of 
the fluctuating nuclear magnetic field correlators. (In the 
limit lo z t c <C 1 we obtain from the above formulas the 
well known result T% = T2, where we took into account 
the fact that the fluctuating nuclear field is isotropic.) 

A more difficult situation arises when lunTq 3> 1, which 
requires a nonperturbative approach. It is this regime 
that we will consider in this paper, i.e. the electron is 
localized in a quantum dot, and the correlation time is 
due to the internal nuclear spin dynamics, i.e., r c = T n 2, 
giving uint c — 10 4 . Thus, no usual treatment and no 
Markov approximation are possible. In particular, the 
perturbative formulas for Ti, T2 (see Eqs. ([!])) are not ap- 
plicable here. Next, we need to address the important is- 
sue of averaging over different nuclear spin configurations 
in a single dot. Without internal nuclear spin dynamics, 
i.e. T„2 — > 00 , no averaging is indicated. However, each 
flip-flop process (due to hyperfine interaction) creates a 
different nuclear spin configuration, and because of the 
spatial variation of the hyperfine coupling constants in- 
side the dot, this leads to a different random value of 
the nuclear field seen by the electron spin and thus to its 
decoherence. Below we will find that this decoherence is 
non-exponential, but still we can indicate the character- 
istic time given by |lOf (A/hN) -1 . Moreover, we shall 
find that T n2 » (A/KN)^ 1 , and thus still no averaging 
over the nuclear configurations is indicated (and dipolar 
interactions will be neglected henceforth). To underline 
the importance of this point, we will contrast below the 
unaveraged correlator with its average. 



II. UNPOLARIZED NUCLEI. 

We consider a single electron in an orbital ground state 
of a quantum dot. Its spin S couples to an external mag- 

I 



netic field B and to nuclear spins {I 1 } via hyperfine con- 
tact interaction, described by the Hamiltonian 

H = gti B S • B + S • hjy , hjv = ^ = gt i B H N , (2) 



where Hjy is the nuclear field. Note that the sum in 
Eq.(H) runs over the entire space. The coupling constant 
with the i-th nucleus, Ai — Avo\^(ri)\ 2 , contains the 
electron envelope wave function ^(rj at the nuclear site 
and vq is the volume of the crystal cell. For sim- 
plicity we consider nuclear spin 1/2. Neglecting dipo- 
lar interactions between the nuclei, we can consider only 
some particular nuclear configuration, described in the 
eigenbasis as \{P Z } >, with P z = ±1/2. |J, [y| Moreover, 
we assume an unpolarized configuration with a typical 
net nuclear magnetic field Aj \\f~N g/is), being much less 
than A/(g[iB) (fully polarized case), for the precise defi- 
nition of TV see below. We study the decay of the electron 
spin from its initial (t = 0) SVeigenstate | f|->. For this 
we evaluate the correlator 



C n {t) =< n I SS z (t)S z I n >, 



(3) 



where <| ... |> means the diagonal matrix element. Here 
5S z {t) = S z {t) - § 2 , and S z (t) = e^p{itn)S z exp(-UH). 
This correlator is proportional to < n\S z (t) — ^(O)!^ >. 
Since at t = the total (electron and nuclear) state 

\n >= I fT,{-^} > is an eigenstate of Hq — S z fiNz 
(with eigenergy e„), we can expand in the perturba- 
tion V = {l/2)(S + h N - + S-h N +) (with H = Hq + 
V). Introducing the usual time evolution operator 
U(t) = f exp(-i J*dtiV(h)), with f being the time- 
ordering operator, we get for the correlator < n \ 
W(t)S z (t)U(t)S z (0) I n >, where the time dependence 
of all operators is due to the H.q Hamiltonian. Then we 
obtain in leading order in V, 



A l 



C n (t) = £l^(cosK fc i) - 1) w 1 V^ [cos(£z + ,/,./)„ + A, /2)l - I] . 

k W nk L £ 2 \ n z)n\ k 



(4) 



where V n k =< n\V\k > is the matrix element between state k =JJ., {••••, 7* = +1/2,...}, and uj nk — e n — e^,, 
initial state n =f|-, {..., if = —1/2,...} and intermediate 
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«z = g^BB z . We have used that \V n k\ 2 — A 2 k < 
n\l/2 - J*|n > /4, and w nfc = e z + (/i z ) n + A fc /2, where 
(h z ) n =< n\h]\f z \n >, and the fact that for the typical nu- 



clear configuration (/i z )„ 



> At. Since N > 1, we 



replace the sums over fc (which run over the entire space) 
by integrals, i.e. J2 k f k = (J d 3 r/v )f(r)+o(l/N). Then 
we have 



Cn(t) 



A 2 



87riV(e z + (^) n )2 



y - 7i(r) cos(e z + (h z ) n )t + 7 2 (T)sin(e z + (^) n )i 



(5) 



r 



where 



/+oo 
-OO 



and 



dzxt{z)F la {Txl{z)), 
(6) 



sin?? cos?7-l sinr? cost; 

= — + — 2 — ; F 2(v) = — • (?) 

Here iV = a z a 2 /vo ^ 1 is the number of nuclei inside 
the dot, and r = At/2irN. Eq.(|) was obtained with the 
use of the following expression for the electron envelope 
wave function: 

\*{p, z)\ 2 = {l/na 2 a z ) eM-p 2 /a 2 )Xo(z)- (8) 

Here a, a z are the dot sizes in the lateral and trans- 
verse (perpendicular to the 2D plane) directions, re- 
spectively, and the transverse wave function Xo( z ) is 

any ana- 

x 2 (o) 



normalized, i.e. f_°° dzXo(z) = 1- For 
lytic function Xo( z ) with expansion Xo( z ) = 
z 2 (Xq)"/2 near its maximum, we have I\,2(j S> 
1) = ±(xg(0)/r 3 / 2 ) v ^VHF[Bin(r X §(0))Tcos(r X 2 (0))]. 
Thus, for this case we obtain the universal power law de- 
cay for times r > 1, i.e. t > (A/N)- 1 : 



C*„(r> 1)~ -7 + 



-3/2 



sia(h n t - 4>), 
(h z ) n +A /2. 



(9) 



Here A is the coupling constant with the nucleus located 
at the center of the dot; and <j> is a phase shift with 
tj> ~ 1. Note that for the typical nuclear configuration the 
quantity A 2 /A^(/i z ) 2 is of order unity, thus, for a weak 
Zeeman field e z < un the part of the electron spin state 
which decays is of the order of the initial value. Hence, 
the same holds for the spin part which survives at r ^> 1 
(i.e. 7 ~ 7 ~ 1/2, for e z < ojn)- We see from Eq.(j|) 
that in the presence of a large Zeeman field, e z 3> loni 
the asymptotic behaviour of C n (t) is not changed, the 
only difference being that the decaying part of the initial 
spin state is small now, i.e. 7 ~ 7 ~ (wjvAz) 2 <C 1. 
For the fully polarized nuclear state Eqs.(^,g) should be 
multiplied by 2, and (h z ) n should be replaced by Aj 2. 

We note that C„(t) in (|J) is quasiperiodic in t, and, 
thus, it will decay only up to the Poincare recurrence time 



Tp. This time can be found from the condition that the 
terms omitted when converting sums to integrals become 
comparable with the integral itself. This will happen at 
r ~ N, giving Tp = 0.1s — Is. 

In next order, V 4 , we face the problem of secular terms 
or "resonances", i.e. the corrections will contain zero de- 
nominators. This gives rise to linearly growing terms 
oc uJNt, even for t -c (A/N)^ 1 . In higher order the de- 
gree of the divergence will increase. This means that the 
decay law we found can, in principle, change after proper 
resummation, because no small expansion parameter ex- 
ists, which, strictly speaking, would justify a perturbative 
approach. Still, the result found in lowest order remains 
qualitatively correct in that it shows that a non-uniform 
hyperfine coupling leads to a non-exponential decay of 
the spin. This conclusion is confirmed by an exactly solv- 
able case to which we turn next. 



III. POLARIZED NUCLEI. EXACTLY 
SOLVABLE CASE. 

In this section we consider the exactly solvable case 
where the initial nuclear spin configuration is fully polar- 
ized (see Fig. 2). We also allow for a magnetic field but 
neglect its effect on the nuclear spins. With the initial 
wave function 



* = I M,T,T,». > 



(10) 



we can construct the exact wave function of the system 
for t > 0, 

*(t) = a(t)¥„ + Yl &WI t, t> ^ t. - >, (11) 



with normalization | ct (i) | 2 +J2k \fik(t)\ 2 = 1, and we as- 
sume that a(t = + ) = l,a(t < 0) = 0. The second 
term in Eq.(|ll|) is an entangled coherent superposition 
of the states with exactly one nuclear spin flipped, and 
thus similar to a magnon excitation. The correlator Cq 
is expressed through a(t) by the formula: 

C (t) = -< ¥ |«&(t)&|*o >= (1 - Ht)\ 2 )/2. (12) 
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Then, inserting ^(t) into the Schrodinger equation we 
obtain 



da(t) 

l : 

dt 
.d0i(t) 



- dt = ( 4 ~ Y m) + Y a{t) + —2~> (13) 

where A = J2k Ak - Laplace transforming (|l3|), a(ui) — 
J °° dt exp(— ujt)a(t), we obtain 



a(uj) 



ia(t = 0+) 



A72- £^1/4(^ + ^/2)' 



(14) 



where iuu — iui — A'/4, A' = A + 2e z , and we have used 
the fact that P k (t = 0) = 0. Now we use the identity 
J2 k A 2 k /A{iu + A k /2) = A/2 - iuj J2 k A k /2{ico + A k /2) 
and replace the sum £ fe iiAj +X k /2 over ^ ne ^JZ-pla 116 by 
an integral. Calculating it using I'J/ (r^) | 2 given in Eq.(||) 
we obtain: 



a{t) 



1 

2nd 



duj i exp[(w - iA'/4)t] 



r iuu + e z + nNicj J dz ln(l - ' 



(15) 

As usual, the integration contour T in Eq. (|15|) is the verti- 
cal line in the complex w-plane so that all singularities of 
the integrand lie to its left (see Fig. 3). These singulari- 
ties are: two branch points (ui = 0,u>o = iAx 2 (0)/2wN = 
iAo/2), and first order poles which lie on the imaginary 
axis (tt> = in). For e 2 > there is one pole, while for 
e z < there are two poles, and for e z — there is one 
first order pole at u>\ w iA/2 + iA J dzxo(z)/4:irN. For 
the contribution from the branch cut between uj = and 
lo = loq we obtain 



e -iA't/4 r , f 

a(t) = - — / dK2z ne tT K < [n 



ttN 



ii 



In 1 1 



K /ttN -2e z /A X 2 Q (0)} 2 + (2tt z ) 2 K 2 } \ 



Xo(z) I 
X 2 o(OW 

(16) 



where r' = txq(0), and zo = zq{k) is defined through 
Xo( z o) = Xo(0) K - We have introduced the dimensionless 
variable k = u>/u>o < 1- In terms of this new variable 
the argument of the log function in Eq. ([l5]) has the form 
1 — Xo( z )/ K Xo(0)- Therefore, for a given k this expression 
changes the sign at Zq(k) which can be found from the 
equation given above. 

Thus, the physical picture can be described as follows. 
At the initial time t = the system has some energy cor- 
responding to the pole and starts to oscillate back and 
forth each time visiting different frequencies within the 
branch cut which corresponds to the flip-flop processes 
with the nuclei located at different sites. Therefore, the 
contribution from the branch cut describes the electron 
spin decoherence. At r of the order of unity (where the 
decay mainly occurs) the decoherence is due to the inter- 
action with the nuclei located at distances of the order of 



the dot radius where the derivative of the coupling con- 
stant is maximal. For longer times, t>1, the asymp- 
totics is determined cither by the interaction with the 
nuclei located far from the dot or near the dot center 
depending on the Zceman field value. 

First we consider the case e 2 = 0, i.e. no magnetic 
field. The asymptotic behavior of Eq.([l6]) for r ^> 1 is 
determined by n <C 1. For example, for Xo( z ) I Xo(0) = 
exp(— z 2 ) we find 



1/ln 3 / 2 ' 



(17) 



This behaviour is not universal and is determined by the 
form of the electron wave function at distances that are 
large compared to the dot size since at large r the deco- 
herence is due to the interaction with the nuclei located 
far from the dot. As a result, a nuclear spin diffusion in- 
duced by the hyperfine interaction occurs in that region. 

Thus, the decay of \a(t)\ starts at r > 1, i.e. at 
t > N/A, as in the unpolarized case (see Fig. 4). Note 
that the magnitude of a is of order of 1/7V (see also 
Rcf.[l6|), thus, the decaying part of the initial spin state 
has this smallness as well, in contrast to the unpolarized 
case above where this part is of order one. The reason for 
this smallness is that for a fully polarized state the gap 
seen by the electron spin through the hyperfine interac- 
tion is A/2, therefore, only a small portion ~ 1/N of the 
opposite (+1/2) state can be admixed. Indeed, in this 
case the change of the energy of the electron subsystem 
is ~ A/N which can be compensated by the energy (order 
of A/N) of a magnon excited in the nuclear subsystem. 

For a large Zeeman field (|e z | 3> ^4) and for r ^ 1, the 
main contribution in ( [l6| ) is given for k — > 1, i.e. by the 
interaction with the nuclei located near the dot center. 
Expanding Xq i z ) for small z (see above), we obtain 
2x§(0)(l - k)/(x§)". 
leJ » A 



Then from Eq.(|16|) we have for 



^^11 -e iT '- iA ' t/4 XM A 2 {l + i) 

a ( T>1 ) = T7=Tt 77=^77 — —Jj2~- ( 18 > 



From this we find then that the correlator ( J12| ) agrees 
with the perturbative result (JsJ) for the fully polarized 
state, i.e. C (t) - C (oo) oc 1/t 3 / 2 . (Note that the 
asymptotic behaviour of the correlator is given by the 
term which is the cross product of the pole contribu- 
tion and the one (Eq. (|l8|) ) from the branch cut.) This 
agreement is to be expected, since for a large Zeeman 
field, the perturbative treatment with a small parame- 
ter A/\e z \ <^ 1 is meaningful (the same is true for any 
model with a small expansion parameter, for example, 
for a system with anisotropy, where the hyperfine con- 
stants in perpendi cular and transverse directions are dif- 
ferent, see section III B|) . However, at zero Zeeman field, 



when the system cannot be treated perturbatively, we 
find Co(t) — Co(oo) oc 1/ ln 3 ^ 2 t, and the agreement with 
(||) breaks down. Nevertheless, the characteristic time 
scale for the onset of the non-exponential decay is the 
same for all cases and given by (A/N) -1 . We have also 
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checked that < S x (t) — S x (0) > has the same behaviour 
as < S z (t) — S z (0) > with the same characteristic time - 
N/A. 

Finally we mention that the observed decay of the elec- 
tron spin can be experimentally studied by local NMR 
measurements IT^l. 



A. Fully polarized nuclei. 2D case 



wave function in different directions. From Eq.([15|) we 
obtain the following singularities: two branch points 
(oj = 0,wo = iA/2wN) and the first order poles which 
lie on the imaginary axis (uj = iv). The position of these 
poles can be found from the equation: 



A/2wNv. (19) 



Here we consider the 2D case when there is no varia- 
tion of the coupling constants in one direction, i.e. we 
use the model representation Xo( z ) = ^(1/2 — \ z \)- This 
case allows us to follow the dependence of the electron 
spin decay law on the spatial variation of the electron 



For positive e z there is only one solution of this equation, 
for negative e z - there are two. Thus, at e z = there is 
one first order pole at ui = ^ [i- eX p(-i/7rA0] • For the 
contribution from the branch cut between to — and 
lu = u>o we obtain: 



J 



a(t) 



exp(-iA't/4) 



dun exp(zr«;) 



[«ln(-l + l//fi) + (/e/7riV) 
I 



(2e z /A)} 2 + tt 2 k, 2 



(20) 



Let us consider first the case of zero Zee- 
man field. Then we have a(t > 0) = 
exp(-^/ 4 )exp(wii)/(7riV[exp(l/7riV) - 1]) + a(t). 
The asymptotic behaviour of Eq.pO) at e z = and for 



r ^> 1 is 1/lnr. In the case of a strong Zeeman field, 
\e z \ ^> A, it follows from Eq.(p0|) that the asymptotics 
at t > 1 is 1/t. Note, that this asymptotics is true in 
the interval Kt« exp(2|e z |/A). Then we obtain: 



J 



< S z (t) > = 



A 2 



2irNe 2 



A 2 



2irNe 2 



Rc { ex P(~ iA ' t / 2 ) ex P( iT ) j 



(21) 



r 



The change of the asymptotics occurs even at \e z \ <C A. If 
1 <C r/lnr <C A/\e z \, then the asymptotics is as before, 
i.e. 1/lnr. For r/ lnr > A/\e z \ > 1 it is 1/(t1ii 2 t). 



B. Anisotropy in the exactly solvable model. 



are different (there is no Zeeman field): Af = A z vq | 
*(r,-) | 2 , A± = A ± v Q I *(r;) | 2 . Then from a system 
of the equations similar to Eq.([l3]) we obtain the fol- 
lowing solution (again for the 2D case, when Xo( z ) = 
0(l/2-|*|)): 



We consider here the model where the hyperfine con- 
stants in the z-direction and in the transverse direction 



J 



duj — 

r iuj 



cxp(~iA z t/4) exp(uit) 



+ (A z /2) - (A 2 ± /2A Z ) + irN(A 2 ± /A 2 z )iu;Hl 



iA z \ 
2ttNu) ' 



(22) 



Note that this equation has the same form as Eq.(15), for the contribution from the branch cut (the decaying 
i.e. an anisotropy plays the role of a Zeeman field. Then part of the initial spin state) we have: 
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exp{-iA z t/4) f 1 



dun exp(iTK) 



ttNv J + + (k/ttTVct) + 1 - l/o-] 2 + 7T 2 K 2 ' 

I 



(23) 



where a = A\/A 2 Z1 and t is defined now as r — A z t/2irN . 
In the case of strong anisotropy, a <C 1, we have at r 3> 1: 



J 



^ 1 2 1 2a 2a exp(-zA^/2) exp(i-r) 
<S z (t)>=--\a\ =--(!-— )-— Re { _ }. 



(24) 



r 



The same result follows for the polarized state from 
Eq.(BI) in the perturbative approach. 



C. Some features of the fully polarized state. 

There are several interesting features which we can ob- 
serve for the fully polarized state. In an external Zee- 
man field, the effective gap seen by the electron spin 
is A'/2 = A/2 + e z . Thus, when e z is made negative 
this gap decreases and even vanishes at \e z \ = A/ 2. 
From Eq.(|l^) we find that the two poles are symmet- 
ric in this case, and the system resonates betwe en the 
two frequencies ui± = ±iA(J Xo(z)dz) 1 / 2 /V&ttN. Note 
that the residual gap is of order A/VN (and not A/N, 
as one might naively expect). Near this Zeeman field we 
have |a(£)| 2 = cos 2 (o; + t) (up to small corrections of order 
l/N), and, as a result, \a\ 2 averaged over time is 1/2, i.e. 
the up and down states of the electron spin are strongly 
coupled via the nuclei (see Fig. 5). In contrast, outside 
this resonance regime the value of \a\ 2 is close to 1 (again 
with small l/N corrections), i.e. < S z (t) >= 1/2 — \a\ 2 
is close to -1/2 at any time. The width of the resonance 
is ~ A/^/N, i.e. small compared to the initial gap A/2. 
We note that this behavior represents periodic (Rabi) os- 
cillations with a single well-defined frequency and is not 
related to decoherence. [The latter is described by the 
branch cut contribution a which remains small (order 
l/N) even near the resonance.] This abrupt change in 
the amplitude of oscillations of < S z (t) > (when chang- 
ing e z in a narrow interval around A/2) can be used for an 
experimental detection of the fully polarized state. Note 
that the weight of the upper pole alone (i.e. that which 
exists at e z = 0) also drops abruptly from a value close 
to 1 to a value much smaller than 1 in the same narrow 
interval, which can be experimentally checked by Fourier 
analysis. 

Another special value of the Zeeman field corresponds 
to the case when the upper pole is close to loq (k = 1) 
the upper edge of the branch cut. This occurs (see 



at the critical value e* 
■dzln|l- X 2 (z)/ X g(0)| < 



= bA/2 < 0, where b = 
— 1 is a non- universal 
number which depends on the dot shape. Since at a finite 
Zeeman field the asymptotics in t is determined by k's 
close to 1, we see from Eq.(fl6|) that for e ~ e* the asymp- 
totics changes abruptly. Indeed, for ((e z — e*)/A) 2 <C 1, 
we find a oc l/y/r, for 1 < r < ((e z - eJ)/A) -2 , and 
a oc 1/t 3 / 2 , for r > ((e z - e*)/A)~ 2 . Thus, when ap- 
proaching the critical Zeeman field e* there is a slow down 
of the asymptotics from 1/t 3 / 2 to 1/t 1 / 2 . It is interesting 
that this slow down is related to a strong modulation of 
the density of states (DOS) of the excitations within the 
continuum band (branch cut) near its edge when e z — * e*. 
In the subspace of none or one nuclear spin flipped (see 
Eq.©), the DOS v becomes 



v(u) = —Im[G (u) + -f-ln-D(w) 

7T du 



(25) 



where u — iu>, Gq(u) = l/(u + A^/2) is the "unper- 
turbed Green's function" , and D (u) is the denominator of 
a(uj), see Eq.(p^). [ pp| The derivation of Eq.(p5|) is given 
in the Appendix. One can then show that for e z — > e* (i.e. 
the upper pole approaches the continuum edge), the DOS 
develops a square root singularity: u(u) oc l/y/uo — u. 
Simultaneously, the weight of the upper pole vanishes 
linearly in t z as e* — e z — > 0. 

Finally, the nuclear state is characterized by /3fc(t), 
which allows for similar evaluation as for a. Here we 
just note that its branch cut part, f3k{t), is nonmono- 
tonic in time, particularly pronounced at t z — > e*: First, 
Pk{t) grows like ^/t, until t reaches ~ 1/(1— dk) 1> an d 
then it decays like l/(y/r(l — au)), with at — Ak/A — > 1. 
Thus, (3k is maximal for A^ close to Aq, i.e. the nuclei 
near the dot-center are affected most by the hyperfine 
interaction with the electron spin. 
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IV. DEPENDENCE OF THE ELECTRON SPIN 
DECOHERENCE ON THE INITIAL NUCLEAR 
STATE. 

So far we have assumed that the initial nuclear state 
had the form of a single tensor product state. Now we 
study the dependence of the electron spin decoherence on 
the initial nuclear state and start with the simple case of 
homogeneous coupling. 



A. Homogeneous coupling. 

In this subsection we consider the case of homogeneous 
coupling, when all the coupling constants are equal, i.e. 
Ak = A/N. It is instructive to start with it, since in 
contrast to naive expectations even in this simple case 
there is some time dependence of < S z (t) > which cannot 
be described by a single frequency On the other hand, 
in this case the expectation value < S z (t) > is a periodic 
function and, for that reason, the electron spin dynamics 
is coherent though the corresponding period can be very 
large for large N (see below). The Hamiltonian is very 
simple, H — (A/N)I-S, where I is the total nuclear spin. 
Then, we can diagonalize H and obtain the result: 

AT 

El= 2N '/ 07 - J i= /+1 /2, 
E 2 = -^^ JorJ 2 = I -1/2. (26) 



In the simplest case of a fully polarized nuclear state and 
opposite initial electron spin polarization (with Wo given 
in Eq.jiq)) we obtain the following result: 



<s z (t) >-- 



2N 



(N + iy 



r(l -cosn N/2 t), (27) 



where fl 



N/2 



A(l + l/N)/2. This result can be eas- 
ily understood, since for the homogeneous coupling I 2 is 
conserved and the initial nuclear state corresponds to the 
maximal value I — N/2. Then from Eq.(26) we obtain 
for the difference of E\ and E% corresponding to / = N/2 
the value which is equal to f2jv/2- Note that the magni- 
tude of the oscillating term is 1/N for N 1, as it has 



been already observed in Sec. Ill 



In the case when initially one nuclear spin is flipped, 
the exact ket has the form 

*i(t) = 5>k(*)U;T,Uk,T,».> + 

k 

+ 2<M*)l1M,Uk,Ui>.-- >, (28) 

k>l 

and the normalization condition is \a.k{t)\ 2 + 
J2k>i \dki(t)\ 2 = 1, with initial condition dki(0) = 0. 
From the solution of the Schrodinger equation we obtain 
for < S z (t) > 



<S z (t) >=- 



1 2(N-2) 

2 + (N - l) 2 



(1 - COsfijV/2-1*) 



2K0)| 2 
N 



(N-2) 
(N-l) 2 

I 



(1 - COsfijV/2-1*) 



2{N - 1 



(N + l) 



2 (1 - COsfl N/2 t) 



where fijv72-i = A{\ — l/N)/2 is the frequency cor- 
responding to the solution E\ — E% of Eq.(|26|) with 
/ = N/2 — 1, and a(0) = afc(0). Now, in contrast to 
the fully polarized case there are oscillating terms with 
two frequencies involved. This is because the initial nu- 
clear state corresponds to I z = N/2 — 1. This value 
of I z can be realized with I = N/2 and I = N/2 - 1. 
Note that for N 3> 1 the amplitude of the oscillating 
term corresponding to / = N/2 is smaller by factor 1/N 
than the amplitude of the oscillating term corresponding 
to I = N/2 — 1. This is a general rule which can be 
checked for an arbitrary initial state with I z = n. This 
state is constructed mostly from the state with I = n 
and the contributions of all other / values are small for 
N > 1. The expectation value < S z (t) >, Eq.(|2|), is 
periodic with period T — AnxN / A. The initial condition 
enters Eq.(^9|) through the quantity a(0) = X)fc afe (0)- 
We can easily check that in the case of the homoge- 
neous coupling there is only a weak dependence on the 



type of the initial nuclear state. For example, the cases 
afe(0) = <5fc n ,a(0) = 1 (which is a single tensor prod- 
uct state) and a randomly correlated (entangled) state 
when the coefficients ctk(0) have random phases, so that 
|a(0)| <C V^V, correspond to the same solution < S z (t) > 
for N 3> 1 , as is easily seen from Eq. (^9|) . 

The solution for the initial nuclear state with a larger 
number of flipped spins can be obtained similarly. For ex- 
ample, when the initial state corresponds to two flipped 
nuclear spins, the solution for < S z (t) > contains oscillat- 
ing terms with three frequencies, Q.n/2, ^/v/2-1 , ^jy/2-2- 
These frequencies are the solutions of Eq.(g6|) for / 
Y 2. X 2 - l,N/2 - 2. Moreover, the initial condi- 
tions (the information about the initial nuclear state) 
enter now < S z (t) > in a more complicated way be- 
sides the quantity |a(0)| 2 which we observed for the case 
with one flipped spin, since it contains also the quan- 
tity J2k l a fc(0)| 2 - Again, depending on the type of the 
correlations between the coefficients «fc(0) we can have 
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different dynamics of the electron spin. However, the de- 
pendence on the type of the initial nuclear state is again 
rather weak. 

B. Inhomogeneous coupling. 

Up to now we have only considered the decoherence 
of the electron spin caused by inhomogeneous coupling 
for a given initial nuclear state which had the form of 
a single tensor product state. Here we consider a more 
general initial nuclear state which can be entangled, i.e. 
contains the coherent superposition of the single tensor 
product states J2r a T\T), where the sum goes over the 
tensor product states \T). This problem was addressed 
numerically in Ref.|l9|. We give here the exact analytical 
solution of a typical problem of this kind which consists 
of a very simple initial nuclear state but still contains 
the relevant physics. In particular, we can examine the 
dependence of the electron spin decoherence (and the cor- 
responding time scales) on the type of initial state. As 
an example of the initial electron and nuclear state which 
we can treat exactly we choose the following one: All but 
one nucleus are polarized in z-direction and the electron 
spin is in the up state. We will see that the results are 
the same for all cases, namely, for this particular nuclear 
state the time scale for the onset of the decay is always 
~ N/A independent of the phases of the coefficients oit- 

The exact ket of the system for t > has the form 

*(t) = X>fe(*)l M,T,U,T,- > +m\ M,T,T,- >, 

k 

(30) 

where the first term at t = represents the initial 
state of the system. The normalization condition is 
Efc KWI 2 + \P(t)\ 2 = 1, and f3(t = 0) = 0. The depen- 
dence on the initial nuclear state comes from a different 
realization of the coefficients a k (0). If, for example, we 



have afc(0) = Skn, then the initial nuclear state consists 
of a single product state. On the other hand, the entan- 
glement of the initial nuclear state can be of a different 
kind, depending on whether the coefficients a k (0) have 
random phases or correlated phases. From the system of 
equations similar to Eq.([l5|) we obtain: 

P{t)=e- aA ' i /|^ex P M)/?(u;), (31) 



„, \ _1 l \ " A k a k (0) 

; [LJ) ~ 2 iw + A/2 - £\ Aj/A{iu> + A/2) + A k /2 ' 

(32) 

Here again A = A k and the integration contour T 
is the same as in Eq. (|l5|) , see Fig. 3. The quantity < 
S z (t) > we are interested in is given by the equation 

< S z (t) >= \ ~ (33) 
As it was already mentioned above, Eq.([32|) contains the 
information about the initial nuclear state through the 
coefficients a k (0)- Note that the denominator of the first 
factor in the solution for (3(uj), see Eq.(|3^), is exactly the 
same as that in a((j), Eq.(n4j), if we put in the latter 
e z = 0. Thus, Eq.(|32|) contains partially the same sin- 
gularities as a(oj). There are, however, some additional 
singularities whose character depends on the properties 
of the coefficients afc(0). 

1. Let us consider the case a k (0) = 5 kn . This means 
that the initial nuclear state is just a simple tensor prod- 
uct state. Then besides the branch cut and the first or- 
der pole (outside the branch cut) which we had before for 
a(u>), there is an additional first order pole lo = iA n /2 
which lies within the branch cut. Considering again the 
continuous limit (i.e. replacing the sum by an integral 
and using Eq.(||)) we obtain: 



Pn(t) 



e -itA/i r duJ 



A 

uit ^-n 



1 



r u" + A n /2) [ 1 + 7rN j dz i n (i _ ^jg) )] 
I 



(34) 



We have evaluated Eq.(|34|) in the 2D case with the model function Xo( z ) introduced in Sec. Ill A and find 
I 



Mr 



_ -itA/4 x 

ttN 



p 



dxe % 



a;(a;o-a;)[7r 2 + ln (1/x-l)] 



JirNr+ir/2 _ MV^O ~ l Y 



z [7r 2 + lr/(l/ar -l)]. 



(35) 



r 



where xq = A n /A < 1, and P means the principle value of the integral. Since /3„(r) is of order 1/N, it follows 
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from Eq.(|33J) that the decaying part of the electron spin 
state has now a smallness ~ 1/N 2 . This is simply due 
to the fact that initially the electron has the same spin 
orientation as all the nuclei except one, and the flip-flop 
process is only allowed with this particular nucleus. How- 
ever, the time scale for the onset of the decoherence is the 
same as before, i.e. it starts at r > 1 and the asymptotic 
dependence of Eq.(|33"| ) at r ^> 1 is 1/lnr, the same as 
before, see Sec. Ill A. 

2. Let us now consider the case when 0^(0) = l/\/~N. 
This corresponds to the entangled initial nuclear state 
where all the terms have the same phases. It is easy to 
see that the singularity of the second factor in Eq. ([32|) (a 
branch cut) coincides with that of the first factor. Thus, 
we obtain that the decay starts at the same time t > 1 as 
it was for the single tensor product state with the same 
decay law, i.e. 1/lnr. 

3. Finally, we consider the case when the phases of the 
coefficients afe(O) are random. It is obvious that while 
the decay law can be different depending on the partic- 
ular choice of the coefficients afe(O) (phases), the char- 
acteristic time scale for the onset of the decay is always 
the same, i.e. N/A. This follows from the fact that the 
singularity of the first factor in Eq.(|32|) (a branch cut) is 
exactly the same as of the second factor. This in turn is 
the consequence of the fact that in the model considered 
here it is only one nuclear spin which is flipped. There- 
fore, the characteristic energy change which determines 
the time scale for the onset of the decay is ~ A/N. 
It is also clear that in a more general case, for example, 
of an unpolarized nuclear state there is a different en- 
ergy scale involved which is related to different values of 
the magnetic fields corresponding to different configura- 
tions which are present in the initial superposition of the 
states. This scale presumably is A/y/N. Then the result 
can be entirely different. It has actually been observed in 
Ref.[l9] that there is a strong dependence of the charac- 
teristic time scale for the onset of the electron spin decay 
on the type of the initial nuclear spin state. 



dtx 

n JO 



<os[y(ti-t 2 )]. (36) 



For r < 1, we get LT^fe cos(A^/2) 
exp[-NC{At/2TrN) 2 }, where C = ir J dz^k{z) / A. 
Thus, the averaged spin correlator C(t), Eq.(|36j), is of 
order — J Q N dx$(x), with $ being the error function. 
Thus, C(t) grows without bound as uj^t for ujpft 3> 1 
(the condition r <C 1 can still be satisfied). Conse- 
quently, the perturbative approach breaks down even in 
leading order in V (we recall that without averaging the 
divergences occur in all higher but not in lowest order). 
To treat this case properly, we need a non-perturbative 
approach. For that purpose we calculate now the 
correlator C(t) exactly by treating the nuclear field 
purely classically, i.e. as a c-number. Then we obtain, 



C„(t) 



h 2 
Ah 2 

N 



(1 — cos h jyt), 



(37) 



where Hn = \fh 2 Nz 



h 2 

jV_L 



h 2 



h 2 



^n± ^ s the nuclear field, with 



The value of corresponds to 
a given nuclear configuration n. To make contact with 
the perturbation procedure we used before in the quan- 
tum case we go to the regime h 2 N , -C h 2 Nz , where /ijv 
can be replaced by Hn z in Eq.Q37j). Then we average 
the resulting expression {h 2 N ^/h^ z )(l — cosh]y z t) over 
a Gaussian distribution for hjsj, i.e. over P(h^) oc 
exp(— 3hj v /2uj 2 j). The result becomes proportional to 
J +o ° dz exp(— z 2 /2)(l — cos(7z))/z 2 cx dx$(x), where 
7 = Thus, we see that we obtain exactly the 

same functional form as before from Eq.(|36"|) with the 
same divergencies in t. This reassures us that the treat- 
ment of the nuclear field as a classical field is not essential. 
On the other hand, the same Gaussian averaging proce- 
dure can now be applied to the non-perturbative form 
Eq.@. Defining C cl (t) = J dh N P(h N )C n (t), we obtain 



V. AVERAGING OVER NUCLEAR 
CONFIGURATIONS. DEPHASING TIME FOR 
AN ENSEMBLE OF DOTS. 



In Sees. |n| and III we have seen that the decay of 
C n (t) occurs in the time interval N/A < t < N 2 /A, 
with N/A ~ 10~ 6 s in GaAs dots. On the other hand, 
the electron spin precesses in the net nuclear field (see 
Eq. (pi)) with the characteristic period (hz)^ 1 ~ lu^ 1 ~ 



1(T W -h lCT 9 s. Thus, lu n l < N/A, and we see that 
the electron spin undergoes many precessions in a given 
nuclear field before decoherence sets in due to the non- 
uniform hyperfine couplings A/-. This behavior changes 
dramatically when we average over nuclear configurations 
JlO| . |2l| ] . For that purpose we consider high temperatures, 
ksT ^> fkjN, and average C n (t) in Eq.([|) over all nuclear 
configurations, i.e. C(t) = J2 n C n {t)/ J2 n - We then find 



C cl (t) 



1 e- 



(38) 



Thus we get rapid (Gaussian) decay of the correlator for 
t ijj/q 1 , giving the dephasing time uj^ 1 — y/N /A. This 
means that < S z (t)S z > saturates at 1/3 of its initial 
value of 1/4. Finally, it seems likely that for the case of 
nuclear quantum spins a non-perturbative treatment of 
the averaged correlator C(t) will lead to a similar rapid 
time decay as found for the classical case in Eq. (|38|). 

VI. CONCLUSION. 

In conclusion, we have studied the spin decoherence of 
an electron confined to a single quantum dot in the pres- 
ence of hyperfine interaction with nuclear spins. The de- 
coherence is due to a non-uniform coupling of the electron 
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spin to nuclei located at different sites. The decoherence 
time is given by HN/A and is of the order of several [is for 
typical GaAs dots. It is shown that in a weak external 
Zeeman field the perturbative treatment of the electron 
spin decoherence is impossible, in particular, we cannot 
use the usual formulas Eq.(|l]). Moreover, the decay of the 
electron spin correlator in time does not have an expo- 
nential character, instead it is given by a power or inverse 
logarithm law. In the case of a strong Zeeman field the 
decay has a universal character oc l/t d ^ 2 , where d is the 
real space dimensionality of the problem. We have also 
solved exactly several model problems which allowed us 
to investigate the dependence of the electron spin decay 
on the initial state of the nuclear system. We have shown 
that there is a strong difference between the decoherence 
time for a single dot, HN/A, and the dephasing time for 
an ensemble of dots, hs/N /A. 



tions: 

and 
where 



PM ~ l u-A>/4 + A l /2 (42) 



ia{t = 0) , i Akpk (t = 0) 

'nA ^—^ i.r7> — A' 



D(Cj) 2D(u))^iu>-A'/4: + Ak/2 



, (43) 



D(u>) =iui + A' /A -iV J^- — 

v ' ' A ^ iuj - A' I A + A k /2 

A k /2 



= iuj + e z + iuj 



and iuo — iu) — A' /A. 

1) Let us start with the case 



E 



iuj + Ak/2 1 



(44) 
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exp(-i«t)|tf >=a(<)|* > +X)j9k(t)|*k >, (45) 

k 

which corresponds to the initial conditions a(t = 0) = 
l,/3fc(i = 0) = 0. Performing the Laplace transform of 
Eq. (B5|) , we obtain 



<* | -\^ >=-ia(Q) = l/D(u>), (46) 

iuj — Ti, 



VIII. APPENDIX 



In this Appendix we derive the density of states (DOS) 
given in Eq.(|2^). The definition for the DOS is 



where we have used Eq.(f43|) with initial conditions a(t = 
0) = 1, and [3 k {t = 0) = 0. 
2) Next we consider the case 

exp{-iHt)\9 k >= a k (t)\^ > >> ( 4 ?) 



f(uj) = —ImTr - 

7T iu-H 



, iu~> = E-ie, e -> 0. (39) 



We write the trace as follows: 

Tr ~^—^ =< ^ol^-^l^o > + E < ^kl-^—^l^k >, 
iuj — H iuj — n , iuj — ti 



where 



(40) 

l*fc>=ltT,T,U,T,->, k = l,...,N, (41) 

The Laplace transform of Eq.(|l3|), a(u>) = 
L dtexp(—Lot)a(t), gives for general initial condi- 



which corresponds to the initial conditions a k (t = 0) = 0, 
and {y£,{t = 0) = Skk'- Again, performing the Laplace 
transform of Eq. ([l7]) , we obtain 

E < >= -*e^)- ( 48 ) 



It follows from Eqs.(^, ^) and the initial conditions 
indicated above that 



AD(iCj-A'/A + A k /2) 2 iuj- A'/A + A k /2' 

^ (49) 
Then from Eqs.(E§), (Ha) we obtain 



E<*^i**>=^E 



y L 



AD ^ (iQ-A'/A+ A k /2) 2 ^ iuj - A' /4 + A k /2 

k k 



(50) 



Collecting all the terms, Eqs . (|4q) , (p0|) 7 we obtain 



where u = ilo. From Eq.([44|) we can easily check that 



v = —Ira 

7T 



k( 1 + E 



D(u>) 



(A k /2f 
{u + A k /2f 



y 1 

Y u + A k/ 2 

(51) 



— InD 

du D 



(A k /2f 
{u + A k /2f 



(52) 
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Thus, from Eq.(51) wc finally obtain Eg. (|25|) of the main text 
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FIG. 1: A single electron spin localized in a quantum dot 
described by a parabolic confinement potential (x,y plane). 
The electron is assumed to be in the orbital ground state 
described by the envelope wave function ^, and interacts with 
the nuclear spins (located at r*;) via hyperfine interaction Ai ~ 
^(ri)! 2 which varies as function of position fi. 
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FIG. 2: The same situation as in Fig.l but with fully polarized 
nuclear spins. 




FIG. 3: The integration contour V in Eq.([L5|), enclosing poles 
and branch cut. 




FIG. 4: Schematic dependence of < S z (t) > on time t for the 
unpolarized tensor product and fully polarized nuclear states. 
The time scale for the onset of the decay ~ N/A is the same 
for both cases. In the fully polarized case the magnitude of 
the effect is 1/N. The period of oscillations is of the order of 
y~N I A for the unpolarized and ~ 1/A for the polarized case. 




FIG. 5: The dependence of |q| 2 averaged over time (< 
|a(i)| 2 >) on the external Zeeman field e z for a fully polarized 
nuclear state. The resonance occurs at \e z \ — A/2, and the 
width of the resonance is ~ A/yN, which is much smaller 
than the initial gap A/2. 



